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ABSTRACT 



04 ' We show that if the intrinsic break energy of Gamma-Ray Bursts (GRBs) is 

determined by the product of more than three random variables the observed 

break energy distribution becomes almost lognormal including the redshift effect 

' because of the central limit theorem. The excess from the lognormal distribution 

at the low break energy is possibly due to the high redshift GRBs. The same 

C\I ■ argument may also apply to the pulse duration, the interval between pulses and 
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so on. 



Among the statistical properties of the observed quantities, lognormal distributions are 
frequently seen in Gamma- Ray Bursts (GRBs). The lognormal distribution may be defined 
as the distribution of a random variable x whose logarithm is normally distributed, 



f(x)dx 



=- exp 

27R7 



(logx-fj,) 2 



dlogx, if x > 0, . . 



0, if x < 0, 



where f(x) is the probability density function for x, and /i and a 2 are the sample mean and 
the variance of logx (e.g., Crow & Shimizu 1988). This distribution is unimodal and posi- 
tively skew. McBreen, Hurley, Long & Metcalfe (1994) pointed out that the total duration 
of the long and short bursts and the time interval between pulses are consistent with the 
lognormal distributions. Li & Fenimore (1996) showed that the pulse fluence and the pulse 
interval distributions within each burst are consistent with lognormal distributions. Nakar 
& Piran (2001) found that the pulse duration also have a lognormal distribution. The break 
energy distribution is also lognormal (Preece et al. 2000; see below). 
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Lloyd, Fryer & Ramirez-Ruiz (2001) suggested that ~ 10% of GRBs might have the 
redshift larger than 6 so that the redshift distribution might be wide. Therefore it is quite 
strange that the observed break energy distribution and the duration distribution are log- 
normal since it does not seem that the observed lognormal distribution reflects the redshift 
distribution of the GRBs. Even if the break energy distribution is lognormal at the source, 
the observed break energy should be smaller than the intrinsic one by a factor of (1 + z) 
while the observed duration should be longer than the intrinsic one by a factor of (1 + z). 
These factors change by order of unity between z = and z ~ 6. 

In this Letter, we will consider a possible origin of the observed lognormal distributions 
in GRBs from the viewpoint of the central limit theorem. 1 

2. BREAK ENERGY DISTRIBUTIONS 

Figure 1 shows the histogram of the break energy taken from the electronic edition 
of Preece et al. (2000). 2 The x 2 test of all data gives the probability of 1.4 x 10 -185 (the 
reduced x 2 is 16.5 with 66 degrees of freedom) that the data was taken from the lognormal 
distribution. Therefore the null hypothesis that the break energy distribution is lognormal 
fails. However, if we exclude the data in the high and low energy ends, the fit becomes 
good as shown in Figure 1. The x 2 test of the data between 70.8 keV and 708 keV gives 
the probability of 0.497 (the reduced x 2 is 0.963 with 17 degrees of freedom) that the data 
was taken from the lognormal distribution, with fi = 2.38 ± 0.004 (Eb — 238 keV) and 
a = 0.240 ± 0.004 (la width is between 137 keV and 413 keV). The improvement of the 
lognormal fit to the break energy distribution excluding the high and low energy ends may 
suggest that the soft and hard bursts originate from a different class of GRBs or emission 
mechanisms. Anyway, hereafter we will assume that the observed break energy distribution 
is lognormal. 

Now let us assume that the intrinsic break energy distribution is lognormal. Next 
we numerically calculate the observed break energy distribution assuming that the redshift 



1 Thc lognormal distribution of the time interval between pulses and the pulse fluence might be reproduced 
by a fine tuning of model parameters (Spada, Panaitescu & Meszaros 2000). 

2 It is not known whether the paucity of the soft and hard bursts is real or not, because harder bursts 
have fewer photons (Cohen, Piran & Narayan 1998; Lloyd & Petrosian 1999; but see Brainerd et al. 1999) 
and there may exist relatively many soft bursts with low luminosities, so called X-ray rich GRBs (or X-ray 
flashes or Fast X-ray transients) (Strohmayer et al. 1998; Heise et al. 2001; Kippcn et al. 2001). Here we 
assume that the selection effect is small. 
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A(l + z) a ~ 1 {\ + z^f^dz, if 0<z<z , 
A(l + z ) a ~ 1 (l + z) h - l d Z) if z <z, 



(2) 



where b < 0. This redshift distribution rises as oc (1 + z) a ~ x to a redshift of z and then 
declines as oc (1 + z) b ~ l , and it is similar to that in Figure 8 of Lloyd, Fryer & Ramirez- Ruiz 
(2001). To mimic Figure 1, we generate 155 bursts with each burst having 35 spectra for each 
realization. We take the mean and the variance of the intrinsic break energy distribution so 
that the mean and the variance of the observed break energy are close to that in Figure 1. 
For each realization, we make a \ 2 fitting to obtain the probability that the data is taken 
from the lognormal distribution. As in Figure 1, we use the data between 70.8 keV and 708 
keV for the \ 2 test- From this simulation, we have found that the average probability can 
reach ~ 0.5 even when the variance of the redshift distribution is comparable to the observed 
one a\ ~ (0.240) 2 , although it is slightly stronger condition to preserve the lognormal form 
than a 2 z < (0.240) 2 . At first glance it is strange that the simulations do not reflect the 
redshift distribution contrary to the argument in Section 1. In the next section we will show 
that it is not strange but natural due to the central limit theorem. 

Figure 2 shows the histogram of the observed break energy for one experimental real- 
ization with (z ,a,b) = (3,3, —3). The x 2 test of the data between 70.8 keV and 708 keV 
gives the probability of 0.128 (the reduced \ 2 is 1-39 with 17 degrees of freedom) that the 
data was taken from the lognormal distribution, with /x = 2.37 ± 0.004 (E b ~ 235 keV) and 
a = 0.256 ± 0.004 (la width is between 130 keV and 423 keV). It is interesting to note the 
excess of soft bursts relative to the lognormal fit as in Figure 1. The average redshifts of 
these soft bursts are relatively high. 



The standard model of the GRB emission is the optically thin synchrotron shock model 
(e.g., Piran 1999). A similar discussion in the following will be applied to the inverse Comp- 
ton model. Let us consider a slow (rapid) shell with a Lorentz factor 7 S (7 r ), a mass m s 
(m r ) and a width l s (l r ). When a separation between two shells is L, the collision takes 
place at a radius R s ~ 2L7 2 . At the collision, the forward and the reverse shock are 
formed. Here we consider the reverse shock propagating into the rapid shell. The discus- 
sion for the forward shock is similar. We assume that a fraction of electrons ( e is acceler- 
ated in the shock to a power law distribution of Lorentz factor 7 e , N{j e )d^ e oc •y~ p d'-y e for 
7e > Imin = [(p — 2)/(p — 1)] (e e u' / ( e n' m e c 2 ) , where n' and u' are the number density and 
the internal energy density in the local frame, respectively, p > 2, and we assume that a 
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fraction e e of the internal energy goes into the electrons. We also assume that a fraction e# 
of the internal energy goes into the magnetic field, B 2 = 87T€bu'. The local frame quantities, 
u' and n', can be calculated using the shock jump conditions (Blandford & Mckee 1976; Sari 
& Piran 1995). We assume that the unshocked shells are cold and the shocked shells are 
extremely hot. If the Lorentz factor of the shocked region is 7, the relative Lorentz factor of 
the unshocked and the shocked region is given by 7 re / ~ (7^/7 + l/lr)/^ — so that 

u ' — {irei — l)n'm p c 2 ~ r y re in'm p c 2 . The number densities of the unshocked and the shocked 
region are given by n' r = m r / Anm p R 2 l r ^ r ~ m r / lQ-Km p L 2 ^ 2 s l r ^ r and n! = (47 rei + 3)n' r ~ 
^1rein' r , respectively. Thus, the characteristic synchrotron energy is given by 



E b = hqeBl1 ^ ~ 260 
m e c(l + z) 



P 



1 2 



p — 1 



X: 



rnl^L-X^XH^l + z)- 1 keV, (3) 



where we assume that the source is at a redshift z. Note that the relative Lorentz factor of 
the unshocked and shocked region 7 re ; depends on the relative Lorentz factor of the rapid 
and slow shell 7 rs ~ (7 r /7s + ls/lr) /2 — 7 r /27 s and the ratio between the number densities 
in these shells / = n' s /n' r = m s l r r ) r /m r l s r ) s (Sari & Piran 1995). For the ultrarelativistic 
shock case 7^ > /, 7re/ = 7 r 1 / 2 / 1/4 /v / 2~ = (m s / r 7^/m r / s 7 s 3 ) 1 / 4 /2. 

Equation (3) shows that the break energy is written in the form of a product of many 
variables. For such a variable made from the product of many variables, the lognormal 
distribution may have a very simple origin, that is, the central limit theorem (Crow & 
Shimizu 1988; Montroll & Shlesinger 1982). Let a variable q be written in the form of a 
product of variables, 



xix 2 



Then, 



log q = log %i + log x 2 H h log x r , 



(4) 



(5) 



When the individual distributions of log Xi satisfy certain weak conditions that include the 
existence of second moments, the central limit theorem is applicable to the variable log q, so 
that the distribution function of logg tends to the normal distribution as n tends to infinity. 

As an example, we numerically generated random variables Xj (i = 1, 2, • • •) whose 
logarithms are uniformly distributed between and 1. Figure 3 shows the histogram of 
the product of these three variables, q = X1X2X3, for 10 4 experimental realizations. The 
distribution of q agrees with the lognormal distribution quite well. It is surprising that the 
X 2 test gives the probability of 0.483 (the reduced x 2 is 1-00 with 278 degrees of freedom) 
that the distribution of q is taken from the lognormal distribution. This example shows 
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that the lognormal distributions may be achieved by a relatively small number of variables 
(Yonetoku & Murakami 2001). Note that when the number of the variables is two, i.e., 
q = x±X2, the probability that the distribution is taken from the lognormal distribution 
was only 1.64 x 10~ 5 , so that a product of only one more variable may make a distribution 
lognormal. 

Therefore the lognormal distribution of the break energy may be a natural result from 
the central limit theorem. We may say, "Astrophysically, not n = 00 but n = 3 gives the 
lognormal distribution !!". The effect of the redshift is just to add one variable in equation 
(3) if the redshift in the observed data is randomly chosen. 

4. PULSE FLUENCE/DURATION/INTERVAL DISTRIBUTIONS 

Let us consider the lognormal distributions in other quantities related to GRBs. When 
the rapid shell catches up the slow one in the internal shock, using the conservation of the 
energy and the momentum, the Lorentz factor of the merged shell 7 m and the internal energy 
Em produced by the collision are given by 7„ t ~ [(m r 'j r + m s ^ s ) / [m r /^ r + m s /7 s )] 1 / 2 and 
Em = m r ^ r + m s ^ s — (m r + m s )^ rni respectively (e.g., Piran 1999). If we assume that a 
fraction e e of the internal energy goes into the electrons and a fraction e w of the energy 
radiated by the electrons is within the gamma-ray band, the observed energy is given by 

E bs = e w e e Eint(l + z)' 1 ~ e w e e m r ^ r (l + z)' 1 . (6) 

Equation (6) shows that the observed energy, which is proportional to the pulse fluence, is 
written in the form of a product of five variables, 3 e w , e e , m r , [y r — r y m +( y rn s /m r )('y s — j m )] ~ 7 r 
and (1 + z)^ 1 . Therefore the lognormal distribution of the pulse fluence may be a natural 
result from the central limit theorem. 

The pulse duration is determined by three time scales: the hydrodynamic time scale, 
the cooling time scale, and the angular spreading time scale (Kobayashi, Piran & Sari 1997; 
Katz 1997; Fenimore, Madras & Sergei 1996). The cooling time scale is usually much shorter 
than the other two time scales in the internal shocks (Sari, Narayan & Piran 1996). The 
hydrodynamic time scale ~ l/c and the angular spreading time scale determine the rise and 
the decay time of the pulse, respectively. Since most observed pulses rise more quickly than 
they decay (Norris et al. 1996), we assume that the pulse duration is mainly determined by 



3 Although the fraction e e of the internal energy that goes into electrons may be a fundamental constant, 
there is still a large dispersion of order unity in e e deduced from the the afterglow observations (Panaitescu 
& Kumar 2001). 
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the angular spreading time ~ R s /2cy, 




(7) 



On the other hand, the interval between pulses At is determined by the separation between 
shells, 



since all shells are moving towards us with almost the speed of light (Kobayashi, Piran & 
Sari 1997; Nakar & Piran 2001). 

Equation (8) shows that the pulse interval At reflects the separation between shells 
L, while equation (7) shows that the pulse duration St is multiplied by one more factor 
(ls/lm) 2 other than (L/c)(l + z). Therefore, if we consider that the distribution of a product 
of variables tends to the lognormal distribution as the number of the multiplied variables 
increases, the distribution of the pulse duration St may be closer to the lognormal distribution 
than that of the pulse interval At. In fact, Nakar & Piran (2001) argued that the pulse 
duration St has the lognormal distribution while the pulse interval At does not, as noticed 
by Li & Fenimore (1996). The pulse interval has an excess of long intervals relative to 
the lognormal distribution. This may suggest the existence of a different distribution, i.e., 
quiescent times (long periods with no activity) (Nakar & Piran 2001; Ramirez- Ruiz & Merloni 
2001). 4 But the central limit theorem may be also responsible for the lognormal distribution 
of the pulse duration St. 



We considered the possible origin of the lognormal distributions in the break energy, the 
pulse fluence and the pulse duration as a result of the central limit theorem. Astrophysically 
the lognormal distribution may be achieved by a product of only a few variables. The effect 
of the redshift is just to add one variable to the product so that the redshift distribution is 
hidden. 

We have no idea about the origin of the lognormal distributions in the pulse interval 5 At 



4 The correlation between the pulse interval St and the pulse duration At is broken when the separation 
between shells is too large to collide each other before the external shock. However, almost all quiescent 

]/ 2 1/3 - 8/3 

times in Nakar & Piran (2001) arc smaller than the limit St < 100 sec E 5 % n 1 ' 7 100 so that shells may 



At ~ (L/c)(l + z), 



(8) 



5. 



DISCUSSIONS 



collide. 



5 Note that the distribution of the pulse interval is lognormal if we exclude the quiescent times (Nakar & 
Piran 2001). 
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and the total duration AT. However the viewing angle may be one factor to be multiplied 
to the pulse interval At. Recently we suggested that the luminosity-lag relation could be 
explained by the variation in the viewing angle 9 V from the axis of the jet (Ioka & Nakamura 
2001; Nakamura 2000). The duration of the pulse from the jet also depends on the viewing 
angle, and according to Figure 2 of Ioka & Nakamura (2001) we have 

Atoc(L/c)(l + ^)(l+ 7 X 2 ), (9) 

when V ~ A6 where 7 is the Lorentz factor of the jet and A8 is the opening half-angle of 
the jet. The multiplied factor (1 +7 2 6^) may be responsible for the lognormal distribution of 
the pulse interval At. Note that the total duration AT is equal to the lifetime of the central 
engine and thus does not depend on 6 V . 
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Fig. 1. — The break energy distribution taken from the electronic edition of Preece et al. 
(2000) is shown. The lognormal fitting of the data between 70.8 keV and 708 keV is shown 
by a solid line. The fitting range between 70.8 keV and 708 keV is shown by dotted lines. 
The x 2 test gives the probability of 0.497 (the reduced x 2 is 0.963 with 17 degrees of freedom) 
that the data was taken from the lognormal distribution, with the mean /x = 2.38 ± 0.004 
(E b ~ 238 keV) and the standard deviation a = 0.240 ± 0.004 (la width is between 137 keV 
and 413 keV). 



500 



z = 3 




10 100 1000 10 4 
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Fig. 2. — One experimental realization of the observed break energy distribution for 
(z , a, b) = (3, 3, —3) in equation (2) is shown. The lognormal fitting of the data between 
70.8 keV and 708 keV is shown by a solid line. The fitting range between 70.8 keV and 708 
keV is shown by dotted lines. The x 2 test gives the probability of 0.128 (the reduced \ 2 is 
1.39 with 17 degrees of freedom) that the data was taken from the lognormal distribution, 
with fj, = 2.37 ± 0.004 (E b ~ 235 keV) and a = 0.256 ± 0.004 (la width is between 130 keV 
and 423 keV). There is an excess of soft bursts relative to the lognormal fit as in Figure 1. 
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Fig. 3. — The distribution of the product of three random variables q = XiX 2 x 3 , each 
distributing uniformly between and 1 in the logarithmic space (that is < logXj < 1), for 
10 4 experimental realizations is shown. The lognormal fitting of the data is shown by a solid 
line. The x 2 test gives the probability of 0.483 (the reduced x 2 is 1-00 with 278 degrees of 
freedom) that the distribution of q is taken from the lognormal distribution. 



